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\x^  Introduction 

Conductive  polymers  such  as  poly(acetylene)  (PA),  poly(pyrrole)  (PP),  and 

poly(paraphenylene)  (PPP)  have  received  much  attention  recently  due  to  their 

potential  applications  and  also  because  of  interest  in  their  fundamental 

properties.  Much  effort  has  been  placed  into  describing  the  mechanism  of  charge 

■/ '•* .  ,  .  transport  through  polymer  films.  Descriptions  of  the  basic  transport  properties 

V  of  conducting  polymers  focussed  initially  on  microscopic  mathematical  transport 

models  where  "electron  hopping*  (1-12]  between  adjacent  oxidized  and  reduced 

sites  in  redox  polymer  films  was  suggested  to  be  the  mechanism  of  charge 

transfer.  Later  it  became  apparent  that  a  consideration  of  concurrent  molecular 

1 

diffusion  of  ions  [13-16]  within  the  polymer  matrix  was  necessary  to  account  for 
observed  polymer  film  behavior.  The  most  popular  treatment  has  involved  coupling 
electron  self -exchange  reactions  within  the  matrix  with  counterion  diffusion  in 
the  polymer  film  [14,  17-26],  This  has  led  to  a  consideration  of  the  measured 
diffusion  coefficient  (Dexp)  as  being  comprised  of  contributions  due  to  ion 
diffusion  (D0)  and  "electron  diffusion"  (Det).  Treated  in  this  manner,  different 
polymer  materials  will  have  different  values  of  electron  transfer  diffusion 
coefficient  and  will  thereby  have  different  conductivities.  To  date  the 
contribution  of  the  charge  transfer  (electron)  diffusion  coefficient  to  overall 
film  transport  behavior  has  been  regarded  as  fundamental  in  describing  the 
transport  characteristics  of  most  conductive  polymers  [27-30].  In  some  cases  the 
contribution  of  "electron  hopping"  to  the  overall  transport  in  various  modified 
electrodes  has  been  determined  to  be  negligible  [31,  32],  and  ion  diffusion 

coefficients  within  polymer  films  have  been  measured  experimentally. 

Most  workers  have  neglected  the  effects  of  migration  on  the  measured  overall 
current.  Descriptions  of  charge  transport  phenomena  have  been  ascribed  primarily 
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to  diffusions!,  effects,  and  have  excluded  other  mass  transport  contributions  to 
the  measured  flux.  Yap  and  Durst  [33,  34]  were  the  first  to  theoretically 
consider  the  effects  of  migration.  Lange  and  Doblhofer  [35,  36]  have 
experimentally  verified  migrational  contributions  in  some  conducting  polymers 
using  a  rotating  ring-disc  electrode,  while  Elliot  [31]  and  Qsteryoung  [32]  have 
observed  migrational  effects  in  some  polymer  electrodes  voltammetrically . 

In  this  paper  a  more  detailed  transport  model  will  be  proposed  which  takes 
into  account  the  migrational  and  diffusional  contributions  to  the  overall 
measured  flux  through  the  film.  Yap  and  Durst  [33,  34]  have  proposed  macroscopic 
models  which  include  migrational  effects,  based  on  the  assumption  of  a  constant 
electric  field  in  the  polymer  film.  In  the  model  presented  herein  a  varied 
potential  gradient  exists  across  the  polymer;  this  must  be  determined  by  the 
local  concentrations  of  all  ionic  species  in  the  solution  within  the  film,  and  by 
the  electrical  conductivity  of  the  film.  The  model  is  based  upon  an  approach 
used  initially  in  describing  the  behavior  of  porous  battery  materials  [37-42], 
In  this  modelling  scheme  the  temporal  behavior  of  the  expected  current  response 
due  to  an  applied  potential  step  is  predicted  theoretically  and  computed 
numerically  via  a  Runge-Kutta  method. 

Experimental 

Numerical  digital  simulations  describing  the  electrochemical  behavior  of 
model  polymer  films  were  carried  out  on  an  IBM  PC/AT.  Two  computer  programs 
(available  from  the  authors)  utilizing  a  fourth-order  Runge-Kutta  algorithm 
combined  with  a  shooting  method  and  an  Adams  predictor-corrector  were  implemented 
(for  similar  algorithms  describing  the  solution  of  coupled  differential 
equations,  see  the  appropriate  sections  in  refs.  [43  and  44]).  In  the  first 
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program,  tha  Initial  conditions  of  the  system  are  set  up,  and  the  local 
parameters  in  the  film  determined  (overpotential,  current,  conversion  as  a 
function  of  position  in  the  film)  for  the  first  time  increment  by  a  shooting 
algorithm.  These  values  are  then  used  as  starting  values  in  the  second  program 
which  evaluates  the  local  overpotential-,  current-,  and  conversion-distance 
profiles  in  subsequent  time  increments.  Simulated  chronoamperograms  are 
constructed  by  integrating  the  resultant  local  current -distance  profiles  and 
examining  the  temporal  behavior  of  the  total  (simulated)  current.  Cottrell  and 
chronocoulometric  plots  are  subsequently  constructed. 

The  Theoretical  Model 

The  Basis  of  the  Model 

It  is  assumed  initially  that  the  charge  (and  discharge)  of  the  model  polymer 
film  is  controlled  purely  by  diffusion  and  migration  of  the  supporting 
electrolyte  within  the  polymer  matrix.  Any  electrochemical  process  that  occurs 
in  the  heterogeneous  polymer  medium  will  be  controlled  by  the  local 
electrochemical  potential,  which  in  turn  will  have  a  time  dependency  determined 
by  the  rate  of  mass  transport  of  the  ions  in  the  medium.  The  Galvani  potential 
of  electrons  in  the  polymer  phase,  ^p,  is  constant,  and  movement  of  charge  is 
fast  compared  with  movement  of  charge  in  the  solution  phase.  Sites  in  the 
polymer  are  converted  at  random  from  neutral  to  oxidized  form  (essentially  a 
Langmuir  model),  with  the  process  being  kinetically  controlled  and  driven  by  the 
difference  in  Galvani  potentials  between  the  polymer  and  solution  phases . 

We  shall  begin  with  the  expression  for  the  local  rate  of  reaction  [45,  46': 


1 


4 

t  a(*r*s)F 

Rate  -  kf’c(o)exp'[ - )  - 

RT 

,  (l-a)(*p-*#>F 

kfe  C(  +  )exp[ - ]  (1) 

RT 

where  kf'  is  the  heterogeneous  rate  constant  for  the  forward  reaction,  kb  is 
that  for  the  reverse  reaction,  a  is  the  transfer  coefficient,  is  the  Galvani 
potential  of  electrons  in  the  solution  phase,  F  is  the  Faraday  constant,  R  is  the 
gas  constant,  T  is  the  absolute  temperature,  and  C(o),  C(+)  represent  the 

concentrations  of  neutral  and  oxidized  forms  of  polymer  within  the  film  matrix, 

i 

respectively.  The  conversion  of  neutral  into  oxidized  form  is  expressed  as 

X  -  C(+)/C°(o)  (2) 

where  C°(o)  is  the  concentration  of  neutral  polymer  in  the  completely  unoxidized 
film.  Then 

(1-X)  -  C(o)/C°(o)  .  (3) 

Inclusion  of  C°(o)  in  the  definition  of  the  rate  constants  yields 

a(*p-*s)F  (l-o)(*p-*s)F 

Rate  -  kfxexp*  [ - ]  -  kb(l-x)exp[ - ]  .  (4) 

RT  RT 

Notice  that  kf  -  kfC°(o)  and  kb  -  kbC°(o). 

We  assume  the  value  of  x  •  0.5  as  the  definition  of  a  standard  state  -  re 
the  difference  in  Galvani  potentials  is  (^p-^3)°.  We  have 


5 


I 


*p  -  *s  -  (*p  -  *s>°  +  n 


(5) 


where  the  overpotential  rj  -  -A^s  and  A^s  is  the  shift  in  the  solution  potential 
from  the  standard  value.  Then 


I 


I 

I 


a(V*s)°F  ar>? 

Rate  -  kf*exp[ - ]exp[ - ]  - 

RT  RT 

(l-a)(*p-rfa)°F  (l-a)f,F 

kb(l-*)exp[ - ]exp[ - ]  .  (6) 

RT  RT 


For  equilibrium  charging  at  a  -  0.5  and  -  0  the  overall  reaction  rate  equals 
zero,  so 


a<V-*s)°F  (l-«)(*p-*s>°F 

kfexp‘[ - ]  -  kbexp[ - ]  -  k°  (7) 


RT 


RT 


and 


arjF  rjF 

Rate-  k°exp'[ - } {** (l-x)exp[ — ])  (8) 

RT  RT 


per  unit  area. 

The  cation  flux  of  the  simple  univalent  electrolyte  is 


3c+  djs 

N+  -  -D+ -  -  c+u+ -  (9) 

3x  3x 
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I 

where  D+  is  Che  diffusion  coefficient  of  the  cation  and  c+  and  u+  are  its 
concentration  and  mobility,  respectively;  x  is  the  distance  coordinate  in  the 
film.  A  similar  expression  can  be  written  to  describe  the  flux  of  anions: 


3c.  3*s 

N.  -  -D. -  +  c.u. - 

dx  dx 


(10) 


Assuming  that  electroneutrality  holds  in  the  solution  phase  (i.e.,  c+  -  c.  -  c), 
we  have  for  the  local  current  density 


3c  d$s 

I  -  F(N+-N_)  -  -F(D+-D.)  —  -  Fc(u++u.) - 

dx  dx 


(ID 


i.e.  , 


d<t> < 


(D+-D.)  3c 


3x  Fc(ul(.+u.)  cCu^+u.)  dx 


(12) 


Substituting  (12)  into  (9), 


3c  u+I  u+(D+-D.)  3c 


N+  -  -D+—  + 


(13a) 


3x  F(u++u.)  u++u_  dx 


3c 

—  -  D  — 

F  dx 


(13b) 


where  the  transporc  number  t+  -  u+/(u++u.)  and  the  binary  diffusion  coefficient  D 


-  2D+D_/(D++D_ ) . 

Alternatively,  N+  can  be  expressed  in  terms  of  I  and  d<t>s/dx.  We  may  now 
take  the  divergence  of  (13)  and,  noting  that 


dl  arjF  rjF 

—  -  aF(Rate)  -aFk°exp'[ - ]  (x- (l-x>exp[ —  ]}  ,  (14) 

5x  RT  RT 


where  a  is  the  specific  area  (equal  to  area/volume  of  the  electroactive 
substrate),  we  arrive  at 


dc  d2c  (i-t+)  ai 

- D  — -  + - .  (15) 

at  ax2  f  ax 


Or  we  may  write 


dc 

a2c  o^f 

rjF 

—  - 

D  -  +  (1- t+)ak°exp"  [ - ] 

{*-(l'X)exp[ — ]) 

(16) 

at 

ax2  RT 

RT 

The  numerical  calculation  is  begun  by  assuming  a  quasi-stationary  state 
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As  these  equations  always  contain  c,  x,  and  rj ,  it  will  be  necessary  to  do  a 
stepwise,  iterative  calculation  [47], 

In  the  quasi-stationary  state  (17)  can  be  integrated  once  to 

dc  (l-t+) 

-  I  ,  <18) 

3x  FD 


while  I  is  given  by  (14).  The  boundary  conditions  at  the  surface  are  I  -  0, 
dc/dx.  -  0,  drj/dx  -  0,  x  -  0,  where  x  -  0  is  the  surface  of  the  polymer  film 
adjacent  to  the  electrode.  The  boundary  conditions  at  the  polymer/  solution 

l 

junction  (x  -  i)  are  c  -  c®,  r,  -  ?appiied  -  -  (*s>applied.  where  c"  is  the  bulk 
concentration  and  £  is  the  thickness  of  the  film. 

Having  lumped  together  the  effects  of  reaction  and  transport  into  one 
equation  to  give  dc/dt  we  must  return  to  the  basic  flux  equation  (11)  and  write 
this 


dc  drj 

I  -  -F(D+-D_) —  +  F(u++u.)c  —  .  (19) 

3x  3x 

Diffusion  coefficients  and  mobilities  within  polymer  electrodes  are  oizen 
experimentally  determinable  [31,32];  if  they  are  unknown.it  may  initially  be 
assumed  that  they  are  given  by  the  solution  phase  values.  In  addition,  we  start 
with  the  initial  conditions  x-0,  c-Ac,  0<x<i;  t-0,  where  Ac  is  a  finite 
quantity  and  the  initial  electrolyte  content  of  the  polymer  is  assumed  to  be 
equal  to  the  bulk  concentration. 


The  scheme  for  the  numerical  calculation  Is  as  follows: 
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i.  An  initial  overpotential  value  is  chosen,  and  the  initial  value  of  xi  is 
calculated  from 

rjiF 

Xi  -  (1-Xi)exp[ - J  -  0  .  (20) 

RT 

ii.  We  choose  the  working  overpotential  and  initiate  the  calculation  by 

substituting  (18)  into  (19)  and  carrying  out  a  simultaneous  integration  of  (14) 
and  (19)  to  satisfy  the  boundary  conditions  at  x  -  0.  Observe  that  if  I  -  0  and 

I 

dq/dx  -  0  at  x  -  0,  then  d^/dx  -  0  automatically  by  virtue  of  (18)  (or  (19)).  I 
is  chosen  at  x  -  i  to  ensure  that  we  obtain  the  boundary  value  at  the  polymer 
film/  solution  interface. 

ili.  Having  derived  I  -  f^(x)  we  derive  c  -  f2(x)  and  x  “  f3(x)  and  use  these 

fixed  values  for  the  next  step  of  the  calculation.  In  determining  x  “  f3(x)  it 

is  assumed  that  I  -  f]/,x)  remains  fixed  over  the  time  interval  At.  Also,  in 

calculating  c  -  f2(x)  we  must  assume  a  value  of  c  at  x  -  l  which  remains  constant 

over  the  duration  of  the  experiment. 

iv.  The  calculation  is  then  extended  to  the  non- steady  state.  Having  found  in. 
steps  i.-iii.)  an  initial  overpotential  which  satisfies  the  bout. '.  irv 

conditions  at  x  -  i,  this  value  can  then  be  introduced  into  a  program  .  n 
calculates  the  conversion-distance  profiles  at  various  times  using  a  non-sn  ,iv 
state  approach.  The  local  currents  and  local  overpotentials  at  varying  dis;  -s 
in  the  film  are  calculated  iteratively,  and  temporal  changes  in  these  cu:  : 
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distance  and  overpotential-distance  profiles  are  computed  using  a  fourth  order 
Runge-Kutta  numerical  method  coupled  with  a  shooting  method  and  Adams  predictor- 
corrector  (43,  44,  48], 


Special  Cases 

Consider  first  the  case  where  slow  mass  transfer  occurs  through  the  solution 
phase  in  the  polymer  matrix,  and  charge  mobility  in  the  polymer  phase  is  very 
fast.  Taking  D+  -  D.  and  a  -  0.5,  we  have 


3»j 

I  -  F(u++u.)c  —  ,  (21) 

dx 


i.e.,  the  concentration  gradient  of  electrolyte  in  the  polymer  film  is 
eliminated.  Differentiating  (21)  and  combining  with  (14)  we  arrive  at 


d2fj  k°  ar/F  r/F 

— -  -  -  exp  ' - {*- (l-x)exp  — ]  -  0  (22) 

dx2  (u++u_)c  RT  RT 


The  solution  of  this  differential  equation  follows  the  same  form  as  the  problem 
of  porous  battery  electrodes  which  has  been  treated  previously  [47,  49],  A  first 
step  in  obtaining  a  solution  is  to  write 


d2rj  Id  dr) 

a»jF 

rjF 

— _  - (  ) 2  -  Texp* 

- (X-(I-X)exp 

-i  , 

(23) 

dx2  2  drj  dx 

RT 

RT 

where  T  -  k°/(u++u_)c.  Integrating,  we  have 
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dr>  , 

RT*r 

at,  F 

RT(l-x)r  (l-a)fjF 

( — )2  - 

- exp' 

— 

- exp -  +  constant  . 

(24) 

dx 

oF 

RT 

(l-a)F  RT 

Note  that  at  x  -  0,  dq/dx  -  0,  and  r,  -  rj0  where  r,0  is  determined  via  the 

previously  described  shooting  method.  Therefore 

RT*r  af?0F  RT(1-X)r  (l-a)»>0F 

constant  -  - exp'  -  +  - exp -  ,  (25) 

aF  RT  (l-a)F  RT 


and 


dr, 

RTXr 

<*r,QF 

ar,F 

—  — 

{ - [exp' 

— 

•  exp'  - ]  + 

dx 

aF 

RT 

RT 

RT(l-*)r  (l-a)*j0F  (l-a)r,F  . 

- [exp -  -  exp - ])^  .  (26) 

(l-a)F  RT  RT 


From  (26) , 


x  - 


1  RTxr  ar?0F  aij'F 

{ - (exp*  -  -  exp*  - ]  + 

°F 


RT 


RT 


RT(l-x)r  (l-a)»?oF  (l-a)r,'F 

- [exp -  -  exp - 

(l-a)F  RT  RT 


])'h  dr,' 


(27a) 


•  l 


arf'F  (l-a)»7'F 

{ab  -  a* exp"  -  +  ef  -  e»exp - 


)'S  dr,' 


RT 


RT 


(27b) 


where  a,  b,  e,  and  f  are  evident  from  the  equation. 

The  strategy  of  the  numerical  calculation  is  as  follows.  r,Q  is  chosen  so  ss 
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to  render  »j  -  »»applied  at  x  -  i  (recall  that  ^applied  is  given  by  the  measured 
current  at  the  polymer/solution  interface  for  an  equipotential  polymer  surface) . 
By  integrating  (27)  we  determine  x  -  f  (q)  and  hence  q  -  g(x)  .  From  this  the 
change  in  x  in  steps  of  time  can  be  calculated.  The  current  is  derived  from  (21) 
at  x  -  I  and,  adding  the  assumption  a  -  0.5,  we  have 


-1 


2RT*r  t?0F  «TF 

( - [exp"  -  -  exp"  - )  + 

„  F  2RT  2RT 

'o 

2RT(1-X)r  »?0f  n'F 

- [exp -  -  exp - ])"4  diy' 

F  2RT  2RT  4 


Now  let  us  consider  the  case  of  low  charge  mobility  in  the  polymer  phase 
coupled  with  high  mobility  of  charge  in  the  solution  phase,  with  the  solution  at 
an  equipotential.  Now  -  constant,  so  rj  -  A^p,  and  the  rate  of  reaction  is 
still  given  by  eq.  (8).  Where  before  it  was  possible  to  ignore  charge  flow  in 
the  polymer  phase ,  in  this  case  the  charge  flow  in  the  solution  phase  may  be 
ignored.  The  charge  flux  in  the  polymer  phase  is  expressed  as 


N+  -  -D, 


observe  that  u+  is  likely  to  be  dependent  on  c+.  A  probable  functional 
relationship  is  electron  transfer  between  occupied  and  unoccupied  sites, 
expressed  macroscopically  as 


u+c+  -  u+'(l-x)c0x 
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where  U+'  is  the  cation  mobility  in  the  unconverted  form  of  the  film,  and  c0  is 
the  initial  electrolyte  concentration.  Now 


3c+  dr, 

I  -  F{  -  D+ -  -  U+'  (l*x)Xco — )  •  (31)  I 

ax  3x 


Applying  the  Nemst-Einstein  relation, 


RT 

D+  -  —  u+  , 

F 

I 

we  have 


RT  3c+  dr) 

I  -  -F( —  u+ -  +  U+’ (l-x)xc0 — )  (33a) 

F  3x  3x 


RT  dX  dr, 

-  -FuV  (l-x)xc0{ - +  — )  ,  (33b) 

Fx  dx  dx 


which  must  be  coupled  with 

i 

* 

I 

ai  qtjF  r,F  I 

- Fk°exp*  [x-d-x)exp— ]  (34) 

dx  RT  RT 

I 

% 

The  boundary  conditions  are  I  -  0,  dri/dx  -  0  at  x  -  2;  r,  -  r,f  at  x  -  0,  with  the 
initial  value  of  x  being  determined  by  t?£.  The  numerical  calculation  is  carried 

! 

out  as  before  except  that  (18)  cannot  be  substituted  into  (33);  rather,  3\  3x  ~ 

must  be  kept  fixed  for  each  step  of  the  calculation. 


► 


(32) 
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Consider  Che  case  where  there  is  low  charge  mobility  in  the  polymer  as  well 
as  low  aobllicy  of  anions  in  the  film.  Now  (29)  is  combined  with  (10)  and  it  is 
assumed  that  electroneutrality  still  holds,  i.e.,  c+  -  c.  -  c.  Two  fluxes  (or 
currents  I+  and  I.)  must  now  be  considered: 


1+ 

F 


(35) 


I.  dc  d+a 

- N.  -  -D.—  +  u.c  -  ,  (36) 

F  3x  3x 

i 

and  the  boundary  conditions  are  3N+/3x  -  0  at  x  -  i,  and  FN+  -  IT  at  x  •  0  (where 
lx  is  the  measured  current  in  the  external  circuit);  dtf./dx  -  0  at  x  -  0,  FN.  - 
It  at  x  -  2.  Observe  that  D+  and  u+  may  be  markedly  dependent  upon  c+  (i.e.  ,  *■)  , 
whereas  D_  and  u.  will  not  vary  appreciably  with  c.  The  anion  flux  can  be 
expressed  as 


dx  d<f>  s 

N.  -  -D.c0 —  +  u.c0*  -  (37) 

dx  dx 


whereas  we  must  use  (30)  or  u+  -  u+'  (where  u+  is  in  this  instance  independer, 
concentration)  to  express  u+  and  D+.  With  (30), 


RT  1  dx 

N+  -  -u+'  (l-x)*c0( - + 

F  x  3x 


d*p 

dx 


(38) 
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The  divergences  of  these  fluxes  will  be  equal  since  the  accumulation  of  (+)  form 
of  the  polymer  (making  no  assumption  as  to  the  identity  of  the  positive  charge, 
be  it  in  the  form  of  polarons,  bipolarons,  or  charged  solitons.)  is  equal  to  the 
accumulation  of  anions,  i.e., 

3N+  3N.  o(*p-*s)F  (*p-*s)F 

- -  -  k°exp*  - [*-(l-X)exp - J  (39) 

dx  dx  RT  RT 


We  may  write 


N. 

D.c0  dx 

d*  s 

4,  ■ 

u.c0x 

u.c0x  dx 

dx 

(40a) 


RT  1  ax  d+s 

- +  -  ,  (40b) 

F  x  3x  dx 


N+ 

_  RT  dX 

3^p 

u+’x(l-x)c0 

Fx  3x 

dx 

Adding  (40)  and  (41), 


1  N.  N+  2RT  dx  dr, 

c0x  u.  u+’(l-x)  Fx  dx  dx 


or 


N_  N+ 

—  +  - 

u.  u+'(l-x) 


2RTc0  dx  dr, 

-  -  c0x  — 

F  dx  3x 


(41) 


(42) 


(43) 
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which  can  be  combined  with  (39)  in  the  usual  form: 


31+  31.  at jF  ijF 

-  - Fk°exp‘  - [X-(1-X)«xp — ]  (44) 

3x  3x  RT  RT 


In  order  to  solve  (43)  with  (44)  note  that  in  addition 


1+  +  I-  -  It  (45) 


FN+  +  FN.  -  IT  ,  *  (46) 


I-  (It-i.) 

dx 

dri 

—  +  - 

-  - 2RTc0 —  - 

Fc0x  “  ■ 

(47) 

u.  u+'U-x) 

dx 

dx 

In  carrying  out  the  iterative  calculation  we  must  "shoot"  to  the  final  value  of 
Ij  which  must  be  adjusted  within  the  calculation  at  each  step. 

Consider  the  case  where  local  equilibrium  is  assumed,  i.e., 


X  rj  F 

-  -  exp — 

l-x  RT 


This  gives 


1  dx  F  r/F  d t) 

-  —  -  — eXp —  — 

(1-X)2  dx  RT  RT  dx 


(48) 


(49) 
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or 


d»j  RT  1  1  dij 

dx  F  exp(^F/RT)  (1-x)2  dx 

RT  1  dX 
F  x(l-*)  dx 

Substituting  into  (47) , 

I.  (IT- I-)  dx  RTc0  dx 

—  +  -  -  -2RTc0 —  - - 

u.  u+'(l-x)  dx  (l-x)tdx 


(50a) 


(50b) 


(51) 


We  now  numerically  integrate,  assuming  as  before  that  x  “  fi(x)  is  fixed  within  a 
given  time  step,  and  evaluate  I.  -  f2(x).  This  distribution  is  then  used  to 
determine  the  next  set  of  values  of  x  “  fi(x). 

Dimensionless  Variables 

In  order  to  solve  the  preceding  problems  it  is  convenient  to  cast 
expressions  in  terms  of  dimensionless  variables.  Recall  that 

dl  art  F  rjF 

—  -  aF(Rate)  -  aFk°exp*  - [*- (l-*)exp — J  ,  (14) 

dx  RT  RT 

which  must  be  solved  with 


I 


dc  drj 

-F(D+-D.) —  +  F(u++u.)c — 
dx  dx 


(19) 
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Taking  Into  account  that 


3c  (l-t+) 

- I  ,  (18) 

3x  FD 


we  have  the  following  dimensionless  variables: 


f  -  —  .  (52) 

RT  » 

I 

1  -  -  ,  (53) 

aFk° 

x 

y  -  -  ,  (54) 

l 

c 

C  -  -  .  (55) 

C(o) 

akTt 

t -  (56) 


Cppp(o) 

Then  (14)  becomes 


3i 

-  exp'{af[x-(l-*)exp(f)] }  (57) 

3y 


For  the  simulation  we  have  the  following  parameters: 


i 
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I 


aFk° 


(58) 


D+C(o) 

7  “  - 

ak°i 


(59) 


(note  that  k°  depends  on  C(o),  i.e.,  7  is  independent  of  C(o)),  and 

I 

D. 

A  —  —  (60) 

D+ 


Observe  that  /9  affects  only  the  current  amplitude ,1  since 


I 

I 


■i 

I  dx 
o 


fi 


(61) 


7  and  A  affect  the  shape  of  the  modeled  current  transient  via  Che  relative 
of  diffusion  and  migration,  along  with  the  volume  rate  of  reaction.  Eq. 
becomes 


rates 

(19) 


ac  af 

i  -  *7(1- A) —  +  7(1+a)C—  ,  (62) 

3y  3y 


while  eq.  (18)  becomes 


ac  i 

ay  27 


1 


(63) 
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The  boundary  conditions  become  i  -  0,  3C/3y  -  0,  3f/3y  -  0  at  y  -  0;  C  -  C®,  f  - 
^applied  at  y  —  1. 

Changes  in  concentration  and  conversion  from  one  step  of  the  calculation  to 
the  next  are  estimated  from 


Ac 

! 


1  31 

- At 

F  3x 


where  At  is  the  time  interval.  This  gives 


ak°  a  l 

Ax  - - At 

C(o)  3y 


or 


3i 

Ax  -  —  Ar 

3y 


where  Af  -  ak°At/C(o) .  Substituting  (63)  into  (62), 


1+A  3f 

K - )  -  7(l+A)C  — 

2  3y 


(6A) 


(65) 


(66) 


(67) 


or 


i 


3f 


2yC  — 
3y 


(68) 
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Thus  the  shape  of  the  expected  current  transient  is  a  function  only  of  X  (and  a) . 

Our  task  is  to  solve  the  simultaneous  equations  (68)  and  (57)  .  Starting 
with  C  -  constant^  and  x  m  constant2.  we  iterate  to  obtain  i  -  f]_(y),  deriving  C 
-  f2(y)  and  x  ~  ^3(y)  •  To  obtain  the  new  value  of  x  we  use  (66),  and  the  new 
concentration  at  each  step  is  determined  from 

fy  L 

—  dy  (69) 

o  2y 

Let  us  apply  a  perturbation,  taking  a  simple  case  where  fD  -  0,  *0  ”  0.5, 
and  C  —  CQ .  We  apply  a  perturbation  f  -  Af,  whens  Af  is  a  small  quantity.  Then 
C  -  C0+AC  and  x  “  Xo+^x.  where  AC  and  Ax  are  functions  of  y.  Since  Af  is  small, 
then  so  is  dAf/dy.  Therefore 

3Af  dAf 

i  -  27(C0+AC) -  -  2yC0 -  ,  (70) 

dy  dy 

and 

d  i 

-  -exp -aAf [ (Xo+^X) *  (1*  ( Xo+A* ) expAf ) ] 

dy 

*  -(1-aAf) [ (x0+Ax) -(l-x0-AX) (1+Af) ] 

-  - (1-aAf ) [ 2Ax*0 . 5Af ] 

=•  -  2Ax  +  0 . 5Af  ( 71 ) 

This  can  be  used  as  an  example  to  examine  the  P.unge-Vutta  numerical  method  r. 
the  first  iterative  cycle  Ax  -  0,  so 
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I 


I 


Differentiating 


di 

-  0.5Af 

ay 


3Af  i 

-  >  mm  ■ 

3y  27C0 


(72) 


(73) 


■ 


yields 


a2Af 

1  3i 

0. 5Af 

77  ’ 

27Cq  3y 

27C0 

(74) 


or 


0.5  0.5 

Af  -  Aexp( - )y  +  Bexp'( - )y  (75) 

27C0  27C0 


At  y  -  0,  3Af/3y  -  0,  or  { (AxO . 5)/(27C0)  -  (BxO. 5)/(27C0) }  -  0,  i.e.,  A  -  B.  A 
y  -  1,  Af  -  Afy«i,  so 

Afy.i  -  A{exp(47C0) _1  +  exp' (47C0) * ^ )  (76) 


Therefore 


2Afy_! 

A  “  -  .  (77) 

cosh[ 1/(47Cq) ] 


i 


Thus 
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■ 

I 


cosh[l/(47C0)]y 

Af  -  Afy_i -  (78a) 

cosh[l/(47C0)] 

l+[l/(47C0)]2y2 

-  Afy_! - —  .  (78b) 

l+[ 1/ (47C0) ] 2 

C  is  regarded  as  a  constant  on  the  next  cycle  (cf.  (73)). 

Results  and  Discussion 

Figures  1-3  (a-c)  show  the  simulated  current-,  conversion-,  and 

l 

overpotential-distance  profiles  for  model  polymer  films  of  thickness  i.  as  a 
function  of  varied  7  (with  a  -  0.5  in  all  cases).  Figures  l-3(c-e)  show  the 
resultant  simulated  profiles  as  the  variable  C  is  varied.  It  is  observed  that  as 
7  is  increased,  the  initial  local  current  increases  with  y,  the  conversion  rate 
increases,  and  the  overpotential  gradient  through  the  polymer  layer  decreases. 
Similar  trends  are  observed  for  increases  in  the  variable  C. 

In  order  to  examine  the  simulated  current-time  transients  expected  during 
film  charging,  the  local  current  at  each  point  in  time  is  integrated  (via  digital 
computer).  This  "measured"  current  (dimensionless)  is  then  plotted  as  a  function 
of  dimensionless  time  r.  Figure  4  shows  the  simulated  current- time  profiles  as  a 
function  of  varied  7  and  C.  It  should  be  noted  that  contributions  due  to 
background  charging  current  and  capacitive  current  have  not  been  included  in 
these  plots.  Departure  from  Cottrell-like  behavior  (i.e.,  pure  diffusion  in  a 
semi -  inf inite  medium,  which  would  predict  i  «  t'1*'  is  due  purely  to  the 
contribution  of  migrational  effects  coupled  with  slow  diffusion  processes  within 
the  polymer  matrix  of  finite  thickness. 
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The  effects  of  migration  [33]  are  seen  more  clearly  in  the  simulated  i  vs. 
t'1*  plots  (Figure  5)  and  in  the  theoretically  predicted  chronocoulometry  (Figure 
6),  where  there  is  observed  significant  departure  from  diffusion  in  a  semi¬ 
infinite  medium  [45,  46].  The  current  vs.  t'15  plots  deviate  considerably  from 
linearity,  and  charge  is  not  proportional  to  J t. 

These  simulated  chronoamperometric  results  differ  from  the  simulations  of 
Yap  and  Durst  [33,  34],  where  it  was  assumed  in  their  calculations  that  a 
constant  electric  field  existed  across  the  width  of  the  polymer  layer.  In  the 
above  simulations,  the  existence  of  a  varying  potential  gradient  across  the  width 
of  the  film  is  included  in  the  model,  thereby  accounting  for  the  resulting 

t 

differences  from  the  simulations  of  Yap  and  Durst  [33]. 

Conclusion 

In  this  paper  a  macroscopic  mathematical  model  has  been  proposed  which 
describes  the  mass  transport  processes  through  conductive  polymer-modified 
electrodes.  The  contributions  of  diffusion  and  migration  have  been  included  in 
the  model;  current-distance,  conversion- dis tance ,  and  overpotential-distance 
profiles  as  well  as  the  overall  simulated  current- time  transients  have  been 
derived  as  a  function  of  the  system  parameters.^-  The  proposed  model  will  be  rhe 
basis  of  other  applications  to  polymer  electrode  systems. 

Observations  of  deviations  from  diffusion- limited  behavior  are  common  in 
conductive  polymer  systems  [31-36,  50-53],  and  models  which  describe  charge 
transport  in  terms  of  diffusion  processes  alone  are  not  sufficient  for  explaining 

^■A  recent  paper  dealing  with  the  transient  response  of  polymer  c  r.ed 
electrodes  under  conditions  of  diffusion  and  migration  of  the  redox  ions  has  ■ .en 
published  since  we  completed  this  manuscript.  These  authors  have  us-  :  ,n 
implicit  finite  difference  technique.  (Rainer  Lange  and  Karl  DoblhofVr  J. 
Electroanal.  Chem.  221  (1987)  13). 


25 


transport  in  many  conducting  polymers.  Hence,  inclusion  of  migrational  and  other 
contributions  into  a  transport  model  {54,  55]  is  often  necessary  in  order  to 
develop  a  more  realistic  view  of  transport  and  kinetic  processes  in  polymer 
electrode  systems. 


I 


I 
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GLOSSARY  OF  SYMBOLS 


a 
a 
0 
c 

Co 
C+ 
c. 
c°° 

c 

C(o) 


Specific  area,  cn'^  (area/volume) 

Transfer  coefficient 
Dimensionless  parameter  (2/aFk°) 

Electrolyte  concentration,  mols  cm'^ 

Initial  electrolyte  concentration,  mols  cm'^ 
Electrolyte  cation  concentration,  mols  cm*^ 
Electrolyte  anion  concentration,  mols  cm*^ 

Bulk  electrolyte  concentration,  mols  cm~^ 
Dimensionless  concentration,  mols  cm*^ 

. .  Concentration  of  neutral  polymer,  mols  cm'^ 

G(+)  Concentration  of  oxidized  polymer,  mols  cm"^ 

C*(o)  Initial  concentration  of  neutral  polymer,  mols  cm'-* 
Conversion  of  polymer  (C(+)/C*(o)) 

Initial  conversion  of  polymer 
Initial  conversion  value 
Binary  diffusion  coefficient,  cm^  s*I 
Electrolyte  cation  diffusion  coefficient,  cm^  s'l 
Electrolyte  anion  diffusion  coefficient,  cm^  <s'I 
Change  in  electrolyte  concentration,  mols  cm*^ 

Change  in  dimensionless  concentration 
Change  in  conversion 
Change  in  dimensionless  overpotential 
Shift  in  Galvani  potential  of  polymer  phase,  V 
Shift  in  Galvani  potential  of  solution  phase,  V 
Time  interval,  s 

Dimensionless  time  interval  (ak”At/C(o) ) 
Overpotential,  V 
Final  overpotential  at  x-0,  V 
Initial  overpotential,  V 
Overpotential  at  x-0,  V 
Dimensionless  overpotential  (»?F/RT) 

Initial  dimensionless  overpotential 
Faraday  constant,  96,484.6  C  eq'I 

Diffuso-kinetic  dimensionless  parameter  (D+C(o)/alk° ) 
Dimensionless  parameter  (k*/(u++u. )c) 

Dimensionless  current  (I/aFi> 

Local  current  density,  A  cm"* 


X 

Xi 

Xo 

D 

D+ 

D. 

Ac 

AC 

Ax 

Af 

A^p 

A*s 

At 

Ar 

n 

If 

H 

f 

h 

F 

7 

r 

i 

I 

1+ 

I. 

It 

kb 

kf 

k* 

kb 

kf 

i 

X 

N+ 

N 

^s 


Current  density  due  to  cation  flux,  A  cm  z 
Current  density  due  to  anion  flux,  A  cm*^ 

Total  current  density,  A  cm'^ 

Concentration  dependent  rate  constant  reverse  rex,  mol  cm^-s 

Concentration  dependent  rate  constant  forward  rex,  mol  cm^-s 

Standard  heterogeneous  rate  constant,  cm  s'l 

Rate  constant  for  reverse  reaction,  cm  s'l 

Rate  constant  for  forward  reaction,  cm  s'l 

Thickness  of  polymer  layer,  cm 

Dimensionless  parameter  (D+/D.) 

Cation  flux,  cm^  s'l 
Anion  flux,  cm^  s'l 

Galvani  potential  of  solution  phase,  V 


I 


J  (kb'C*(o) 
1  (kf ' C° (o) 
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Galvani  potential  of  polymer  phase,  V 
R  Molar  gas  constant,  8.3144  J  mol'^-K'^ 

t  Tine ,  a 

t+  Cation  transport  number  (u+/(u++u.)) 

t.  Anion  transport  number  (u./Cu^+u.)) 

T  Temperature,  K 

r  Dimensionless  time  (ak*t/C<o)> 

u*.  Electrolyte  cation  mobility,  cm/s 

u.  Electrolyte  anion  mobility,  cm/s 

x  Distance  from  electrode,  cm 

y  Dimensionless  distance  parameter  (x/i) 


Figure  1.  Current  vs.  distance  profiles  for  varied  7  and  C;  a-0.5.  (a) 

7-1.0.C-0.5;  (b)  7-IO,  C-0.5;  (c)  7-100,  C-0.5;  (d)  7-100,  C-0.3;  (e) 

7-100, C-0. 15.  Time  increases  from  top  to  bottom  for  r-0  to  60  in  equal 

increments . 

Figure  2.  Conversion  vs.  distance  profiles  for  varied  7  and  C;  a-0.5.  (a)  7-I.0, 

C-0.5;  (b)  7-IO,  C-0.5;  (c)  7-IOO,  C-0.5;  (d)  7-IOO,  C-0.3; (e)  7-IOO, 
C-0. 15.  Time  increases  from  bottom  to  top  for  r-0  to  60  in  equal 
increments . 

Figure  3.  Overpotential  vs.  distance  profiles  for  varied  7  and  C;  a-0.5.  (a) 

7-1.0,  C-0.5;  (b)  7-IO,  C-0.5;  (c)  7-100,  C-0.5;  (d)  7-IOO,  C-0.3; (e) 
7-100,  C-0. 15.  Time  increases  from  bottom  to  top  for  r-0  to  60  in 

equal  increments. 

Figure  4.  Simulated  chronoamperograms  for  varied  7  and  C.  (a)  7  varied  at 
C-0. 15;  (b)  C  varied  at  7-50.  * 

Figure  5.  Simulated  current  vs.  l/^t  plots  for  varied  7  and  C.  (a)  7  varied  at 
C-0. 15;  (b)  C  varied  at  t-50. 

Figure  6.  Simulated  chronocoulometric  plots  for  varied  7  and  C.  (a)  7  varied  at 
C-0. 15;  (b)  C  varied  at  7-50. 
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